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In an earlier paper we have analytically determined the photon regions and the shadows of 
black holes of the Plebahski class of metrics which are also known as the Kerr-Newman- 
NUT—(anti-)deSitter metrics. These metrics are characterized by six parameters: mass, 
spin, electric and magnetic charge, gravitomagnetic NUT charge, and the cosmological 
constant. Here we extend this analysis to the Plebahski—Demiahski class of metrics which 
contains, in addition to these six parameters, the so-called acceleration parameter. All 
these metrics are axially symmetric and stationary type D solutions to the Einstein- 
Maxwell equations with a cosmological constant. We derive analytical formulas for the 
photon regions (i.e., for the regions that contain spherical lightlike geodesics) and for 
the boundary curve of the shadow as it is seen by an observer at Boyer-Lindquist coor¬ 
dinates (ro,'ho) in the domain of outer communication. Whereas all relevant formulas 
are derived for the whole Plebahski—Demiahski class, we concentrate on the accelerated 
Kerr metric (i.e., only mass, spin and acceleration parameter are different from zero) 
when discussing the influence of the acceleration parameter on the photon region and 
on the shadow in terms of pictures. The accelerated Kerr metric is also known as the 
rotating C-metric. We discuss how our analytical formulas can be used for calculating 
the horizontal and vertical angular diameters of the shadow and we estimate these values 
for the black holes at the center of our Galaxy and at the center of M87. 
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1. Introduction 

Basically, the shadow of a black hole is the region on the observer’s sky that is left 
dark if the light sources are anywhere in the universe but not between the observer 
and the black hole. For a mathematical description, it is convenient to consider 
light rays that are sent from the observer’s position into the past. Some of them 
are deflected by the black hole and then go out to bigger radii again; because they 
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can reach one of the light sources, we assign brightness to their initial directions. 
Others stay close to the radial line and go to the horizon without meeting one of 
the light sources; to their initial directions we assign darkness. The resulting dark 
region on the sky is called the shadow of the black hole. The boundary of the shadow 
is determined by light rays that spiral towards a lightlike geodesic which stays on 
a sphere. The region outside of the black hole filled with these spherical lightlike 
geodesics is called the exterior photon region. Recently, the shadow of a black has 
attracted even Hollywood’s attention. In the movie Interstellar, which was released 
in 2014, it was shown how nearby observers would see the shadow of an almost 
extremal rotating (Kerr) black hole with an accretion disk.l^ 

In the near future astronomers actually expect to observe the shadow of a black 
hole. Currently there are two cooperating projects—the US-led Event Horizon Tele¬ 
scope projeciP and the European BlackHoleCam project—who try to image the 
shadow of galactic black holes. It has been predicted since several years that there 
are supermassive black holes at the centers of most—if not all—galaxies. There is 
strong evidence for a black hole at the center of our own Galaxy, associated with 
the radio source Sagittarius A* (Sgr A*): Infrared observations of orbits of stars 
near the centeilll0 demonstrate that there has to be a heavy object with a mass of 
approximately 4.3 million Solar masses within a small volume. The most convincing 
candidate for such an object is a black hole. Another good candidate for a supermas¬ 
sive black hole is the object at the center of M87 (m > 10® Solar masses). Because 
of the large distance, the diameter of the shadow of galactic black holes will be tiny. 
Even for the nearest candidate Sgr A* (8.3 kpc awaj00), the predicted diameter of 
the shadow is only about 50microarcseconds, see Sec.[^ Although tiny, such a diam¬ 
eter should be resolvable with very large baseline interferometry (VLBI) 0111 Numer¬ 
ical studies taking scattering into account suggest that the shadow can be observed 
only at sub-millimeter wavelengths, see Falcke, Melia and Agol.^^ The resolution of 
interferometric measurements in this wavelength regime will be further improved if 
the 10-meter space-based radio-telescope Millimetror^ goes into operation, proba¬ 
bly in the mid-2020s. With this Russian satellite the Earth-based telescope network 
is upgraded with an extra-long Space-Earth baseline of 1.5 million km. 

If the shadow of a black hole will be observed, its shape will give important 
information on the parameters of the black hole. Therefore, it is necessary to pro¬ 
vide the observers with theoretical calculations of the shape of the shadow for 
different black-hole models. In an earlier paper,lil we have given an analytic for¬ 
mula for the shadow of black holes of the Plebahski class, which are also known 
as Kerr-Newman-NUT-(anti-)de Sitter space-times. Black holes of this class are 
characterized by mass, spin, electric or magnetic charge, a NUT-parameter and a 
cosmological constant. In the present paper we extend this analysis to the bigger 
class of Plebahski-Demianski space-times, which include in addition a so-called ac¬ 
celeration parameter. The Plebahski-Demiahski space-times are axially symmetric 
stationary solutions to the Einstein-Maxwell equations with a cosmological constant 
and they are of Petrov type D. Before determining the boundary curve of the shadow 
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in these space-times, we have to study the photon region. We develop the relevant 
formulas for the whole class of Plebahski-Demiahski space-times. However, when 
studying the effect of the acceleration parameter onto the photon region and onto 
the shadow in terms of pictures, we restrict to black holes for which only the mass, 
the spin and the acceleration parameter are different from zero. This is because 
the effect of the other parameters - electric and magnetic charge, NUT-parameter 
and cosmological constant - has been studied already in our earlier paper and the 
presence of the acceleration does not change this significantly. The metric of these 
black hole space-times characterized by mass, spin and acceleration alone is known 
as the rotating C-metric or the accelerated Kerr metric. 

If only the mass and the acceleration parameter are different from zero, we have 
the so-called C-metric which describes a space-time with boost-rotation symme¬ 
try. This solution to the vacuum Einstein field equation was found by Levi-Civita 
( 1919 jiSI and Weyl (1919)in]llll The name C-metric refers to the classification in 
the review of Ehlers and Kundt (1962).E^ The rotating version of the C-metric was 
considered by Hong and TecP^ while a detailed discussion of accelerated space-times 
in general can be found in the book by Griffiths and Podolsky.E^ 

Commonly the C-metric is given in the form introduced by Hong and TecJ^^ 

C ^ ^ 

with cubic functions G = (I — + 2amx) and F = —(1 — — 2amy). The 

metric depends on two parameters, the mass m and the acceleration parameter a. 
The domain covered by the coordinates (r, x, j/, ft) actually contains two black holes 
accelerating away from each other with a conical singularity (a “strut”) on the axis 
of rotational symmetryFor our purposes, Boyer-Lindquist coordinates are 
more suitable, see Eq. ([^ below, which cover only one of the two black holes. 

There are several earlier papers on the shadows of black holes. Here we just 
mention some important works, for a more comprehensive list we refer to the in¬ 
troduction of our earlier paper.^^^ The first analytic calculations of the shadow of 
a black hole were done by Syng^^ for the Schwarzschild metric (Synge used the 
word “escape cone” for the complement of the shadow) and by BardeerPi^ for the 
Kerr metric. In the Schwarzschild case the photon region reduces to the “photon 
sphere” at r = 3m and the shadow is circular. In the Kerr case, the photon region 
is spatially three-dimensional and the shadow is non-circular. The deviation of the 
shadow from a circle could be used as a measure for the spin.l^ With ray-tracing 
algorithms it is possible to include effects of matter on the light rays like emission 
regions and scattering.l^SlMI Here we do not take such effects into account but re¬ 
strict ourselves to the purely geometric construction of the shadow based on the 
assumption that light rays are lightlike geodesics and that there are no light sources 
between us and the black hole. It is our goal to derive an analytical formula for the 
shape of the shadow from which, in future work, the parameters of the black hole 
could be extracted with analytical means. 
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After a discussion of relevant properties of the Plebahski-Demiahski space-times 
(Sec. we determine the photon region (Sec. which is essential for calculating 
the boundary of the shadow of the black hole for an observer at Boyer-Lindquist co¬ 
ordinates (tqt'&o) (Sec.[^. We derive all relevant formulas for the whole Plebahski- 
Demiahski class. However, when illustrating the results with pictures of the photon 
region and of the shadow in Secs. andwe restrict to the accelerated Kerr metric. 
In Sec. we discuss how our analytical formulas can be used for calculating the 
angular diameters of the shadow. We use these equations for estimating the angular 
diameters of the shadows of Sgr A* and M87. 


2. The Plebahski Demiahski metrics 

We consider the general Plebahski-Demiahski class of stationary, axially symmet¬ 
ric type D solutions of the Einstein-Maxwell equations with a cosmological con¬ 
stant. As a matter of fact, these solutions were first found by DebeveJ^ in 1971 
but are better known in the form of Plebahski and DemiahskPi^ from 1976. The 
Plebahski-Demiahski metrics are discussed in detail by Griffiths and Podolskj^iS 
and by Stephani et al.^^ It is common to use rescaled units, i.e. units so that the 
speed of light and the gravitational constant are normalized (c = 1, G = 1). In 
Boyer-Lindquist coordinates this metric is then given by, see Ref. [T^ p. 311 

da;'" = ^ dr^ -f ^ dr?^) -b ^ ((S -h ax)^A^ sin^ 'd - dtp^ 

2 / 

-b — (^ArX - a(S -b ax) A^ sin^ d 
with the abbreviations 

n = l- g(^-bacosi?)r, 

Tj = + (£ + acos-d)^, 

X = asin^ d — 2^(cosi? -b G), 

The coefficients of the polynomials and A^ are 

03 = 2^am - 4a£{^{k + /3) + f), 

04 = -b /3) -b f), 

6o = fc -b /3, 

hi = —2m, 

b2 = ^+ - (o^ + 3f){^ik + /?) + f), 

^3 = -2S (A - + /3) + t))), 

h = -{^k+^) 


= 1 — 03 cos 19 — 04 cos^ l 9 , 

Ar = 6o + bir + + bsr^ + h^r'^. ^ ^ 



( 5 ) 
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with 


k = 




l + 3^eia^-P) 




oj = \ld?- A- P‘- (6) 


Basically, the coordinates t and r may range over all of K while d and ip are 
standard coordinates on the two-sphere. Note, however, that for some values of the 
black-hole parameters r and 'd have to be restricted, see below. The Plebahski- 
Demiahski space-time depends on seven parameters, to, a, /3, £, a, A and C, which 
are to be interpreted in the following way. to is the mass of the black hole and a is 
its spin. /3 is a parameter that comprises electric and magnetic charge, /? = 9e +9mj 
if non-negative; if j3 is negative, the metric cannot be interpreted as a solution to 
the Einstein-Maxwell equations but metrics of this form occur in some brane-world 
scenarios.^ The NUT parameter £ is to be interpreted as a gravitomagnetic charge. 
The parameter a gives the acceleration of the black hole which is at the center 
of our investigation while A is the cosmological constant. The parameter C, which 
was introduced by Manko and Ruiz,l^ is relevant only if £ ^ 0. In this case there 
is a (conic) singularity somewhere on the z axis and by choosing C appropriately 
this singularity can be distributed symmetrically or asymmetrically on the positive 
and the negative z axis. Note that this parameter C has nothing to do with the 
name “C-metric” for the accelerated Schwarzschild space-time. All the parameters, 
TO, a, £, /3, A, a and C, may take arbitrary real values in principle, albeit not all 
possibilities are physically relevant. 

If a = 0, the Plebaiiski-Demiahski class reduces to the Plebahski space-timed^ 
which are also known as the Kerr-Newman-NUT-(anti-)de Sitter space-times. For 
this more special class of metrics the photon regions and the shadows have been 
discussed in our earlier paper, see Ref. The Plebahski-Demiahski class covers 
many well-known non-accelerated (a = 0) space-times like the Schwarzschild (a = 
/3 = £ = A = 0), Kerr (/3 = £ = A = 0), or Reissner-Nordstrom space-time 
(o = £ = A = 0) as well as the accelerated C-metric (a = /I = £ = A = 0) or their 
rotating version (/? = £ = A = 0) which we will call accelerated Kerr space-time. 

The metric I© becomes singular at the roots of U, E, A^, and sind. Some 
of them are mere coordinate singularities while others are true (curvature) singu¬ 
larities. In the following we briefly discuss the influence of a on these singularities, 
becomes zero if 


^ s/a? 

q;(£-I- acosd)' 

As the metric blows up if U —0, Eq. (0 determines the boundary of the space- 
time, i.e., we have to restrict to the region where U is positive. The allowed region 
is a half-space bounded by a plane (£ = 0), a half-space bounded by one sheet of 
a two-sheeted hyperboloid (£^ < a?), a domain bounded by a cylinder (£^ = cd), 
or a domain bounded by an ellipsoid (£^ > a^), see Fig. For a = 0 there is no 
restriction because U = I. 
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S becomes zero at the ring singularity 

r = 0 and cos'd = —£/a. ( 8 ) 

It exists for P' < a? and is a curvature singularity (if m 7 ^ 0). Outside of this 
singularity the sphere r = 0 is regular, so it is possible to travel through one of 
these two hemispheres (“throats”) from the region r > 0 to the region r < 0 and 
vice versa. If P > a?', there is no ring singularity and the sphere r = 0 is regular 
everywhere. In the limiting case where = a?' the ring singularity degenerates into 
a point on the axis. It becomes a point singularity for (. = a = Q that disconnects 
the space-time into the regions r > 0 and r < 0. The ring singularity is unaffected 
by a. 

Moreover, the metric is singular on the z axis, i.e. if sin'd = 0. If a 7 ^ 0 oder 
f 7 ^ 0, this is a true (conical) singularity on (at least a part of) the rotational axis. In 
the NUT case the singularity depends on the Manko-Ruiz parameter C. For further 
details we refer to the book by Griffiths and Podolsky.li^ 

The real roots of A,, yield coordinate singularities which correspond to the up 
to 4 horizons ri > r 2 > ... of the space-time. If a = 0 and A = 0, then A,, reduces 
to a second-degree polynomial, A^ = — 2mr -\- a? — + £3^ and horizons can be 

found at 

= m ± P — P (9) 

if < a^^y. := m? + P — j3\ then r+(= n) is the outer (event) horizon of the 
black hole and r-{= r 2 ) is the inner horizon. The domain of outer communication 
is the region outside of the outer horizon where A^ > 0. For we would 

find, instead of a black hole, a naked singularity or a regular space-time. Since we 
are interested only in the black hole case, we will not consider this possibility in 
the following. In the accelerated or cosmological scenario (a 7 ^ 0 or A 7 ^ 0) the 
horizons could not in general be specified in a simple form because A^ is then a 
polynomial of degree 4. Depending on the sign of the leading coefficient 64 , which 
is mostly determined by whether > P and by the sign of A, the vector field 
dr is timelike or spacelike for big values of r, see Fig. In the timelike case (left 
column in Fig. [^, all real roots of A^ are in the allowed region with U > 0. Hence, 
the first root ri represents a cosmological horizon and the subsequent root r 2 is the 
black-hole horizon. In this case the domain of outer communication is the region 
between ri and r 2 where A,. > 0 (gray shaded and hatched region in Fig. [^. In the 
case that dr is spacelike for big r (right column in Fig. [^, the first root ri is not 
in the allowed region. Hence, we have at r 2 a cosmological horizon, at the event 
horizon of the black hole, and in between the outer domain of communication. The 
horizons can be easily determined if /3 = £ = A = 0. Then k = of and to = a, hence 

Ar = (r^ — 2mr + a^)(l — a^r^) ( 10 ) 

and we find the usual (Kerr) horizons at r = r± given by ([^ with i = 0 and /? = 0, 
and the additional horizons at r = ±-. Of course, we must have |a| < —. 

Oi ■ II r+ 
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Fig. 1. A schematic illustration of the graph of Ar (upper row) and a polar plot of the region 
n > 0 (lower row) given by 0 - Depending on the sign of the leading coefficient 64 , see A, 
goes to +00 of —00 for big radii r; the space-times (with A > 0) belonging to the plots in the 
left column are dominated by the Kerr property (£^ < a^) and in the right column by the NUT 
property (a^ < £^). The space-time is restricted to that region where Q > 0 which is marked here 
with a light gray shading. Geometrically, the boundary of this region is one sheet of a two-sheeted 
hyperboloid (left) or an ellipsoid (right). As in the NUT dominated case (right) the root ri of Ar 
is not inside the allowed ellipsoid Q > 0, it could be no event horizon. Thus, the event horizon of 
the black hole is instead at rs (right) or at r 2 (left). The gray-shaded and hatched region marks 
the outer domain of communication (Ar > 0 ) where we will place the observers. 


In general, at the roots of would be coordinate singularities, too; these would 
indicate further horizons where the vector field would change the causal character 
from spacelike to timelike, just as the vector held dr does at the roots of Ar- 
However, since these horizons would lie on cones d = constant instead of on spheres 
r = constant such a situation would be hardly of any physical relevance. Therefore, 
we exclude it by limiting the parameters of the black hole appropriately. As = 0 
implies 


cosd± = 


—03 ± + 404 

204 


( 11 ) 


0 is guaranteed for all real if the radicand in is negative. Therefore, we 
agree to choose the black-hole parameters such that 0 > O 3 -I- 404 . If /3 = £ = A = 0, 
this condition can be simplihed. Then k = and lo = a, hence 


A-s = 1 — 2am cos -I- a^o^ cos^ 


and A^ 7 ^ 0 is assured if 


a < 


2m 

r_ _ 

7 ^ — 


( 12 ) 


m — \/m?' — d- 


if O = 0 , 
if O 7 ^ 0 


(13) 
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There are some other interesting regions around a black hole characterized by 
the change of the causal character of the Killing vector fields dt and d^p. 

In that region where dt becomes spacelike, i.e. gtt > 0, no observer can move on 
a t-line. Thus, any observer in this region has to rotate (in ip direction). This region 
with gtt > 0 is known as the ergosphere or the ergoregior^ An ergoregion only exists 
if a 7 ^ 0. Note that at the horizons, i.e. at the roots of A^, the metric coefficient 
gtt = (Ar — sin^ i?) is positive. Hence, the horizons are always contained 
within the ergoregion. For a 7 ^ 0 or A 7 ^ 0 there are “cosmological” horizons in 
addition to the black-hole horizons; then the ergoregion consists of several connected 
components. The boundary of (a component of) the ergoregion is always tangential 
to the horizon on the rotational axis, i. e. at r? = 0 , tt. 

If a 7 ^ 0 or 7 ^ 0, there are regions where the Killing field d^, becomes timelike, 
g^pip = 0. This indicates causality violation, because the tp-lines are closed timelike 
curves. For £ = 0 the region where g^^ = 0 is completely contained in the domain 
where r < 0 and, thus, hidden behind the horizon for an observer in the domain 
of outer communication. In the case I' 7 ^ 0, however, there is a causality violating 
region in the domain of outer communication around the axial singularity. 

In the following, we will only consider the black-hole case where we have at 
least one positive root of A^. Observers will be placed in the domain of outer 
communication. 


3. Photon Regions 

In the Plebahski class of space-times, i.e., for a = 0, the geodesic equation is com¬ 
pletely integrable; in addition to the obvious constants of motion, there is a fourth 
constant of motion, known as the Carter constant, which is associated with a second- 
rank Killing tensor. If a 7 ^ 0, instead of this Killing tensor we only have a confor¬ 
mal Killing tensor. This is sufhcient to assure complete integrability for lightlike 
geodesics. The four constants of motion are the Lagrangian 

C = \gp,vX^x'', £ = 0 for light, (14) 


the energy E and the z-component of the angular momentum 


^ dC . . 

E ■ = -^ = 

at 


dC 

9‘^tt-, ( 1 ^) 


and the Carter constant see Ref. 36, which is now associated only with a confor¬ 
mal Killing tensor. The Carter constant may be viewed as the separation constant 
for the r and the ^ motion of lightlike geodesics. The four constants of motion allow 


^Some authors call only the intersection of the region where gtt > 0 with the domain of outer 
communication the ergoregion. This is that part of the region gtt > 0 which an outside observer 
would be able to see. 
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US to write the lightlike geodesic equation in separated hrst-order form, 

X(L^ - Ex) (E + ax)((S + ax)E-aL^) 

sin^ d Ar ’ 

Lz - Ex ^ a{i^ + ax)E - aLz) 

A.g sin^ d Ar ’ 

=:eW, 

Sin 17 

^ ^ - ArK =: R{t). 



E . 


E 

11^^ = 


)V = 

' E 


112 



(16a) 

(16b) 

(16c) 

(16d) 


In order to derive an equation for the shadow of accelerated black holes, we proceed 
in the same way as for the Plebahski space-times. As a first step, we have to de¬ 
termine the spherical lightlike geodesics, i.e., lightlike geodesics staying on a sphere 
r = constant. The region filled by these spherical geodesics is called the photon 
region K,. Mathematically, spherical orbits are characterized by f = 0 and r = 0 
which requires by ( |16d ) that R{r) = 0 and R'{r) = 0. Thus 

{(T. + ax) - aLsy Ar{{J: + ax) - aL e) 

Ke = -- , Ke = --, (17) 

where A), is the derivative of A^ with respect to r and Ke, Le are abbreviations 

(18) 
Hz 

After solving ([T7| for the constants of motion 




Le = —. 
E 


Ke = 


16r^Ar 

jKyp 


oLe = (E -I- ax) - 


ArAr 


(19) 


we can substitute these expressions into (16cI. As the left-hand side of (16c) is 
non-negative, we find an inequality that determines the photon region 


A: (^ArAr — EA),)^ < Ida^r^A^A^ sin^ i?. 


( 20 ) 


Of course, the equality sign determines the boundary of the photon region. 

Just as in the non-accelerated space-time£ll2Zl for every point (rp, dp) of /C there 
is a lightlike geodesic through {vp, dp) that stays on the sphere r = rp. The d motion 
is an oscillation bounded by the boundary of JC while the ip motion given by (16b) 
might be rather complicated. 

The stability of these spherical geodesic with respect to radial perturbations is 
determined by the sign of i?"; a spherical geodesic at r = rp is unstable if R''{rp) > 0, 
and stable if R"{rp) < 0. From (16dl we get with (19) 


R"{r),„ 


8E^ 


A? = 2rA,AL 


^A'J - 2r^ArA: 


( 21 ) 
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A non-rotating black hole (a = 0) is surrounded by a photon sphere^ rather than 
by a photon region, since the inequality ( |20| ) defining 1C reduces to an equality 

( 22 ) 


ArAr = (r^ +£2 )a;. 


The best known example is the photon sphere at r = 3m in the Schwarzschild 
space-time. 

Because of the rotational symmetry it is convenient to plot a meridional section 
through space-time for illustrating the regions around a black hole. The resulting 
pictures, which are {r,'d) polar diagrams where § is measured from the positive z- 
axis, are shown in Fig. Each figure contains the photon region /C, where unstable 
and stable spherical light rays are distinguished according to (21), the horizons r± 


of the black hole given as boundaries of the region where < 0, the ergoregion, 
the causality violating region, and the ring singularity. 

The dashed circle marks the throats at the sphere r = 0. For viewing the whole 
range of the space-time, we use two different scales for the radial coordinate: In 
the inner region r < 0 (inside the sphere r = 0) the radial coordinate is plotted as 
mexp (r/m); this is continuously extended with r -|- m in the outer region r > 0 
(outside the sphere r = 0). By not plotting just the exponential of the Boyer- 
Lindquist coordinate r, as suggested by O’Neill,!^ we avoid a strong deformation of 
the outer part. 

While our formulas apply to black holes of the entire Plebaiiski-Demianski class, 
in the pictures we restrict to accelerated Kerr space-times (/3 = £ = A = 0) because 
we want to focus on the effect of the acceleration. Fig. comprises various images 
of photon regions /C for different spin and acceleration parameters, where the spin 
is varied in the rows and the acceleration in the columns. For the spin we choose 
fractions of the value for an extremal black hole a = Aomax with As 
and Oniax = rn, cf. ([^ and (10). Although one would expect only very small accel¬ 
eration parameters in reality, we choose relatively big values (a € {O, 
for a better illustration of the effects. For each of the three values of the acceleration 
the figure is compared with the ordinary Kerr case. 

In the ordinary Kerr space-time (a = 0), see the left half of the images in the 
columns of Fig.[^ we see two photon regions: one with unstable orbits in the exterior 
region of the black hole at r > r_|_ and one in the interior region at r < r_ which 
contains unstable orbits as well as stable ones. For spinning black holes the exterior 
photon region develops a crescent-shaped cross-section which grows with increasing 
spin a. The inner photon region consists of two parts divided by the ring singularity. 
Note that also in the rotating case there are circular photon orbits, namely at that 
five points on the boundary of /C which are tangent to a sphere r = constant: there 
are three circular photon orbits in the equatorial plane—two at the boundary of the 
exterior photon region and one at the boundary of the interior photon region—and 
two more off the equatorial plane at the boundary of the interior photon region 
where r < 0. Furthermore, we find the ergoregion containing the horizons of the 
black hole, and in the interior adjacent to the ring singularity a causality violating 
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region with < 0 
boundary (A^ = 0): horizons 
unstable spherical light-rays in /C 
stable spherical light-rays in K, 


^ causality violation < 0) 
^ ergoregion {gtt > 0) 
throats at r = 0 
• ring singularity 



Fig. 2. Photon regions in accelerated Kerr space-time for spins a = Aamax, where amax = In 
each column the plot for the unaccelerated Kerr space-time (left) is compared to the plot for an 
accelerated Kerr space-time (right). The specific acceleration parameters are listed in the bottom 
row. A legend for the plotted regions is given at the top. 


region. If > 77i^/2 the ergoregion intersects the exterior photon region. All of 
these regions are symmetric with respect to the equatorial plane. 
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The plots for the accelerated Kerr space-times look similar to the non-accelerated 
ones but there are two significant differences. Firstly, a non-zero acceleration param¬ 
eter gives rise to additional horizons, similarly to a cosmological constant. Secondly, 
the plots are no longer symmetric with respect to the equatorial plane which is sim¬ 
ilar to the NUT case. The additional outer horizon, a cosmological one, is best seen 
in the illustration for the highest acceleration a = In principle, such a horizon 
also appears in all other plots but most or even all of it is located outside of the 
shown clipping. The asymmetry with respect to the equatorial plane is best seen 
for a = ^. With the exception of the causality violating region, the entire picture 
looks as if pushed into the negative z direction, i.e., into the direction against the 
direction of the acceleration. For a better view. Fig. shows bigger versions of the 
two plots shown in the fourth row of the third column in Fig. 

As one would expect, the photon region, the ergoregion and the causality vio¬ 
lating region depend on the signs of a and a. While the photon region is reflected 
at the equatorial plane if the sign of a is changed, the ergoregion and the causality 
violating region are reflected if the sign of aa is changed. The effects of (3, i and A 
on the photon regions have been discussed in our earlier paper, see Ref. We do 
not repeat this here because there are no new qualitative aspects if a is present. 



Fig. 3. Bigger versions of both plots shown in the fourth row of the third column in Fig. 
Illustrated are the photon regions of spinning Kerr black holes (a = gm) where the left plot belongs 
to the ordinary space-time (a = 0) and the right plot to an accelerated space-time (a = 


4. Shadows of Black Holes 

If one looks into the direction of a black hole then there is a region on the sky which 
stays dark, provided that there are no light sources between the observer and the 
black hole. This dark region is called the shadow of the black hole. To determine 
the shape of the shadow we consider light rays which are sent into the past from the 
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position (ro,i!^o) of a fixecj^ observer in the domain of outer communication. Then 
we can distinguish between two types of lightlike geodesics: Those where the radial 
coordinate increases after possibly passing through a minimum and those where the 
radial coordinate decreases until reaching the horizon at r = r_|_. If we assume that 
there are light sources distributed in the universe, but not between the observer and 
the black hole, geodesics of the first kind could reach a light source; so we assign 
brightness to the initial direction of such a light ray. Correspondingly, we assign 
darkness to the initial directions of light rays of the second kind, i.e., these initial 
directions determine the shadow of the black hole. The boundary of the shadow 
corresponds to light rays on the borderline between the two kinds. These are light 
rays that spiral asymptotically towards one of the unstable spherical light orbits in 
the exterior photon region /C. Hence, the essential information for determining the 
shadow of a black hole is in the surrounding photon region. One may even say that 
the shadow is an image of the photon region (but not of the event horizon). 

For deriving an analytical formula for the boundary curve of the shadow we 
proceed, again, as in the case without acceleration. First, we choose an orthonormal 
tetrad, cf. page 307 in Ref. 15 for our fixed observer at (tq, do) 

(E + ax)dt + ad^ 

(ro,’do) 


Co = H 


\/E Af 


62 = —H 


+ xdt) 


-/EA^jsind 




ei — 



eg = -H 


(ro,‘do) 



(23) 


(ro,‘do) 


Since our observer is in the domain of outer communication, is positive. E is 
positive everywhere (except at the ring singularity which is not part of the space- 
time and, moreover, away from the domain of outer communication) and A^ is 
positive by assumption. This guarantees real coefficients in Eqs. ( |2^ . It is easy to 
check that the Ci are orthonormal. As usual, the timelike vector Cq is interpreted as 
the four-velocity of our observer. By our choice of the tetrad, eg ± eg are tangential 
to the principal null congruences of our metric; eg points into the spatial direction 
towards the center of the black hole. So we have chosen the four-velocity of our 
observer adapted to the symmetries of the space-time in the sense that eg is in 
the intersection of the t-tp-plane and the plane spanned by the two principal null 
directions. 

For any light ray A(s) = (r(s), d(s), :/?(s), t(s)), the tangent vector at the position 
of the observer can be written in two different ways, using either the coordinate basis 
or the tetrad introduced above. 


X = fdr +'dd,} + + idt, (24) 

A = tT( — eo -f sindcos^ei -|- sin d sin'(/'eg + cos deg). (25) 


^Because of the symmetry, it is enough to specify the r and •& coordinate to define a fixed position 
in space-time. 
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observer at (ro,^?o) 



Fig. 4. Eq. ( |25[ l defines celestial coordinates 0 and ijj for the light rays at the observer’s position, 
as illustrated in the left figure. With this choice, 0 = 0 is the direction towards the black hole. 
Every direction of a light ray represented by a point {0,'ip) on the celestial sphere (black ball) is 
visualized by its stereographic projection (red ball) on a plane. The dotted (red) circles are the 
celestial equator 9 = Trf 2 and its image under stereographic projection. 


The second equation defines the celestial coordinates 9 and ip, see Fig. For the 
scalar factor a we obtain with (TT^ 


= 5(A,eo) = 


aLz - (E + ax)E 




(26) 


iro,‘So) 


We substitute and r from (16b I and (16d) into (24), and we insert the expressions 
of Ci from (23) into (25). Then comparing the coefficients of d^p and dr in the 
resulting two equations yields 

T :=sin6» = 


X^ArKE 


P :=smtj} = 


J.2 £2 _ uLe 

Le + a cos^ -9 + 2£ cos 9 


(27a) 


(27b) 


'd—9o 


\/A^Ke sind 

where 

LE = LE-a + 2iC. (28) 

If a light ray asymptotically approaches a spherical lightlike geodesic at a radius 
Tp in the photon region, it must have the same constants of motion as this limiting 
spherical geodesic. By (19) and (28), this implies that the constants of motion of 
light rays that correspond to boundary points of the shadow are given by 


KEirp) = 


16r^A,- 


(a;)= 


aZrir,) = 


(29) 
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Here the range of Vp is determined by the intersection of the exterior photon region 
(20) with the cone § = § 0 ^ cf. Ref. Thus, for a rotating black hole Vp ranges over 
an interval at whose boundary points (|20|) holds for = -do with equality. If we 


insert (29) into (27), we get the boundary curve of the shadow on the observer’s sky 


parametrized with Vp. In the case a = 0 the photon region degenerates into a photon 
sphere r = Vp. This unique value defines a unique KE{rp) but does not restrict 


Le- Calculating the corresponding 9 from (27) gives the radius of the shadow which 


is circular in this case. We may use Lb as a parameter for the boundary curve. 


where Le varies between the extremal values given by (16c) for 0(i?o) = 0. 

Comparison with Ref. shows that the formula (|20[) for the photon region as 
well as the formulas (27 28 29) for the boundary curve of the shadow are identical 


with those of the non-accelerated case. However, the metric functions (|^ have now 
a more general meaning because they include the acceleration parameter. 

Several properties of the shadow are preserved, even with the acceleration pa¬ 


rameter added. A non-rotating black hole still has a circular shadow since (27a) 


depends on the unique KE{rp) but not on Le, so 9 = constant in this case. As the 
acceleration parameter breaks the spherical symmetry, this is a non-trivial result. 
Furthermore, the shadow is still independent of the Manko-Ruiz parameter C which 
is relevant only in the case £ ^ 0. As in the non-accelerated case, the shadow is al¬ 
ways symmetric with respect to a horizontal axis, because {tp, 9) and (tt — ip, 9) are 
determined by the same constants of motion Ke and Le- Again, this is a non-trivial 
result because it is not implied by an underlying symmetry unless .^ = 0, a = 0 and 
■do = 7r/2. 

It is to be emphasized that we have calculated the shape of the shadow for an 
observer with a particular four-velocity, adapted to the principal null directions of 
the space-time. For an observer in a different state of motion, the shadow is distorted 


by aberration. These aberration effects have been discussed in detail in Ref. 39 As 
the aberration formula maps circles onto circles, the statement that a non-rotating 
black hole produces a circular shadow is true for an observer in any state of motion. 

Figures and comprise images of shadows for different space-times seen by 
an observer at tq = 3.8m with varying inclination -do- As explained in Fig. 
we map the shadow onto a plane by stereographic projection. Standard Cartesian 
coordinates in that plane of projection are given by 


x{p) 

y{p) 


= —2 tan {\9{p)) 


sin'0(p) 

cosip{p) 


(30) 


In Fig. we show the shadow for accelerated Kerr space-times where we have 
chosen the same values for a and a as in Fig. Here, the observer is fixed at 
Boyer-Lindquist coordinates ro = 3.8m and ■do = 7r/2 (in the domain of outer 
communication). The different values of a are encoded into different shadings. 

Also with acceleration, the shape of the shadow is largely determined by the 
spin a of the black hole. Hence, the shadow becomes more and more asymmetric 
with respect to a vertical axis with increasing spin a where the asymmetry results 
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Fig. 5. Shadows of accelerated Kerr black holes (omax = rn) seen by an observer at rp = 3.8m 
and = 7r/2 for different spin values. The magnitude of the acceleration is color-coded where 
the specific values of q o are listed below the plots. The dashed (red) circle marks the projection 
of the celestial equator, cf. Fig. 

from the “dragging effect” of the rotation on the light rays. The shadow is reflected 
at a vertical axis if the sign of a, i. e. the spin direction, is changed. One might 
have expected a similar effect with respect to a horizontal axis if the sign of a is 
changed. However, this is not true. As the shadow stays symmetric with respect to 
a horizontal axis even if a yf 0, the shadow is independent of the direction of the 
acceleration, i. e. of the sign of a. The acceleration has an effect on the size of the 
shadow, as is visible with the naked eye. This, however, has little relevance in view 
of observations because the size also scales with rQ and a comparison of the radius 
coordinates in different space-times has no direct operational meaning. 
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= f = f ^0 = 1 1^0 = 0 


Fig. 6. Shadow of a black hole in accelerated Kerr space-time {a = a = m = amax) for an 
observer at vq = 3.8m with different inclination angles As in Fig.^ the dashed (red) circle 
indicates the celestial equator. 


With the plots in Fig. we investigate the influence of the observer’s inclination 
on the shadow of an extremal Kerr black hole (a = Omax = w) with acceleration 
a = As in Fig. the observer ist fixed at tq = 3.8to. Clearly, for —>• 0 the 

shadow becomes circular. We have already emphasized the remarkable fact that the 
shadow is always symmetric with respect to the horizontal axis. 


5. Angular Diameters of the Shadow of Black Holes 


From the analytical formulas (271 and (29) for the boundary curve of the shadow 
we can deduce expressions for the horizontal and vertical angular diameters of the 
shadow. These correspond to the dashed lines in Fig. Owing to the symmetry, 
the angular diameters 5h and are determined by three angular radii Qh ^, , and 

Qv as indicated in Fig. 

5h = Qhi+Qh 2 , sin = sin sin 0,,, = (31) 

6t, = 2gy, sinp„ = costjjv sin9)v = \/l — P^(rt,)T(r„), (32) 


where T and P have the same meaning as in ( [27| ). 

In the following we restrict to the Kerr space-time with an observer in the 
equatorial plane, 'do = f • Even in this case, a formula for the angular diameters 



Fig. 7. Angular radii of the shadow of a black 
hole. Owing to the symmetry with respect to 
a horizontal axis, the two angular diameters 
(dashed lines) of the shadow are given by three 
angular radii: two horizontal radii q^. and one 
vertical radius The angular diameters are 
calculated as 5^ = and 5v = 2^^,, 

respectively. 
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of the shadow was not known before, as far as we know. In the general case, the 
angular diameters can be calculated analogously; it is true that then the radius 
values Thi and are zeros of a polynomial of higher than fourth order, so they 
cannot be determined in closed form. In terms of these radii, however, one gets 
analytical formulas for the angular diameters also in the general case. The horizontal 
angular radii gh^ are characterized by iphi = ±f, so we must solve the equation 


1 = sim'0(r/i) = P^{rh) which in the Kerr case simplifies to (use Eq. (27b) with 


Eq. (29)) 


rh{rh 

— 3m) ^ = 4ma^ 

(33) 


rhi = 2m + 2m cos(C/3), 

(34a) 


rh 2 = 2m — mcos(C/3) — •\/3msin(()/3), 

(34b) 


rh 3 = 2m — mcos(C/3) -I- •\/3msin(()/3), 

(34c) 


where ^ = arg((2a^TO — m^) — i{2amy/'m? — a^)). Here we have to choose the solu¬ 
tions and rh 2 which are the radii of the two circular photon orbits in the exterior 


photon region. Evaluating PT for rti and yields by (31) the horizontal angular 
diameter 6h of the shadow. 

The vertical angular radius corresponds to those boundary points where the 
tangent is horizontal. By (32) we have /(r„) := sin^ p.;, = — p2(r„))T^(r„), so 

the tangent is horizontal if^^(r„) = 0. This yields 

0 = (1 - P^)T' - PP'T\r^ (35) 

_ y^A(ro) Ty (r„(2a^ -|- Tq) - imry) {a?m - r«(3m^ - 3mr„ -|- r^)) 


A{ry{ry{2a^ + ’’o + ’’5) “ + 3’’^)) 

where we have to choose the unique solution inside the exterior photon region 


(36) 


SmrQ 


2a2 

With this value r„ we get an analytic expression of the vertical angular radius 

& = (1 - P^)T%. =- 27mVS(.^ + ro(ro-2m)) - 


(37) 


(38) 


Tq -I- QoPtq -I- 3a2(4a2 — Q'w?)r‘Q -\- 8a® 

For a = 0, we recover from (38) Synge’s formulcP^for a Schwarzschild black hole. 


2 27m? (ro — 2m) 

sm g = -o-. 


(39) 


Since the shadow of a non-rotating black hole is always circular, the horizontal 
angular radii gti are also given by (39) in this casej^Note that for all values 0 < 
< m^ (38) gives the same value as (39), 27m?I tq, if m is negligibly small in 


'^For a = 0 one finds = arg(—m®) = —tt and ^ = "im. Then T^{3m) reproduces ||39|l. 
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comparison to tq- This means that for observers far away from the black hole the 
vertical diameter of the shadow is independent of a. 

In the extremal Kerr space-time, a = m, the circular photon orbits are at rh-^ = 


Am and rh 2 = w since C, = arg(m^) = 0. Together with (37) this results in the 
following formulas for the angular radii 

64m^(ro — 


sm Qh^ = 


sm Qh 2 = 


(r^ + 1 


: Q -T 8m2)2 


sm Qy = 


27m^rQ 


(ro -I- mY{r‘^ + Sm^)' 


(40) 


(ro -I- m)2 


Finally, we use (39) and (40) to determine the angular diameters given by (31) 


and (32) for the shadow of the black hole in the center of our Galaxy near Sgr A* 


and of that in M87. The resulting values are given in Table together with the 
corresponding values for the mass M (in multiples of the Solar mass Mq) and the 
distance ro of the black holes. We use two sets of parameters for M87 because the 
mass estimation based on the modeling of stellar dynamics yields a mass twice as 
big as the estimation based on gas dynamical measurements, compare Refs. [40] - |^ 
The horizontal angular diameter for the extremal rotating black holes is al¬ 
ways about 13% smaller than for the Schwarzschild case while the vertical angular 
diameters Sy coincide in all cases. We have already observed that the latter is a 
consequence of the fact that ro is large in comparison to m. It turns out that the 
shadow of the black hole in M87 is not much smaller than that of the black hole at 
the center of our Galaxy; the bigger distance of M87 is almost compensated by its 
bigger mass. 


Table 1. Horizontal and vertical angular diameter Sy of the shadow for Sgr A* 
and M87 for a non-rotating Schwarzschild model (a = 0) or an extremal rotating Kerr 
model (a = m) of their black holes 



Sgr A* 

M87 

Sh 

M87 

Sh 

a = 0 

a = m 

53.1 lias 
46.0 lias 

53.1 lias 
53.1 yas 

37.8 lias 37.8 pas 

32.8 lias 37.8 nas 

20.1 lias 20.1 lias 
17.4 lias 20.1 lias 

m=MG 

M = 4.31 X 10® Mq, 
ro = 8.33 kpdmi 

M = 6.2 X 10® Mq, 
ro = 16.68 

M = 3.5 X 10® Mq, 
ro - 17.9 MpJl£lEll 


6. Conclusions and Outlook 

We have seen that knowing the photon region surrounding a black hole is essential 
for calculating the shadow. For both the photon region and the boundary curve 
of the shadow we have found analytical formulas in the general type D class of 
PlebahskEDemiaiiski space-times. Since these space-times are not in general asymp- 
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totically fiat and possess additional cosmological horizons, it is not possible to re¬ 
strict to observers at infinity as it was done in many other articles on shadows of 
black holes. We have placed our observer at any Boyer-Lindquist coordinates in 
the outer domain of communication instead, and we have calculated the shadow 
for the case that the four-velocity of the observer is adapted to the symmetries of 
the space-time in the sense that it lies in the intersection of the t-(p-plane with the 
plane spanned by the principal null directions. Interestingly, for such an observer 
the shadow is always symmetric with respect to a horizontal axis, independently of 
an acceleration a 7 ^ 0, an inclination 7 ^ 7 r /2 of the observer, or a gravitomagnetic 
NUT charge £ 7 ^ 0. The boundary curve of the shadow depends on the space-time 
parameters m, a, £, /3, a, and A, as well as on the observer’s position {ro,'do)- For 
an observer whose four-velocity is not adapted to the symmetries of the space-time, 
the boundary curve of the shadow also depends on the 3 components of the spatial 
velocity with respect to our standard observer. 

Although the acceleration parameter does not destroy the symmetry of the 
shadow with respect to a horizontal axis, it does have such an effect on the pho¬ 
ton region, the ergosphere and the causality violating region. The photon region is 
reflected at the equatorial plane if the sign of a is changed whereas the ergoregion 
and the causality violating region are reflected at the equatorial plane if the sign of 
aa is changed. 

Our estimates of the angular diameters for the shadows of the black holes in the 
centers of our Galaxy and of M87 show that the shadows are roughly of the same 
size. Hence the planned observations may provide us with shadow images not only 
of the black hole in our Galaxy but also of that in M87. If the current attempts 
of observing the shadow are successful, this will give us a chance to deduce the 
parameters of the black hole from the boundary curve of the shadow. Our analytical 
formula combined with a Fourier analysis should be a promising tool for achieving 
this goal. We are planning to investigate this in a follow-up article. 
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